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ABSTRACT 

We present archival Rossi X-ray Timing Explorer (RXTE) and simultaneous Advanced 
Satellite for Cosmology and Astrophysics {ASCA) data of the eclipsing low mass X- 
ray binary (LMXB) X 1822—371. Our spectral analysis shows that a variety of simple 
models can fit the spectra relatively well. Of these models, we explore two in detail 
through phase resolved fits. These two models represent the case of a very optically 
thick and a very optically thin corona. While systematic residuals remain at high 
energies, the overall spectral shape is well-approximated. The same two basic models 
are fit to the X-ray light curve, which shows sinusoidal modulations interpreted as 
absorption by an opaque disc rim of varying height. The geometry we infer from these 
fits is consistent with previous studies: the disc rim reaches out to the tidal truncation 
radius, while the radius of the corona (approximated as spherical) is very close to 
the circularization radius. Timing analysis of the RXTE data shows a time lag from 
hard to soft consistent with the coronal size inferred from the fits. Neither the spectra 
nor the light curve fits allow us to rule out either model, leaving a key ingredient of 
the X 1822—371 puzzle unsolved. Furthermore, while previous studies were consistent 
with the central object being a 1.4 M© neutron star, which has been adopted as the 
best guess scenario for this system, our light curve fits show that a white dwarf or 
black hole primary can work just as well. Based on previously published estimates of 
the orbital evolution of X 1822—371, however, we suggest that this system contains 
either a neutron star or a low mass 2.5 Mq) black hole and is in a transitional 
state of duration shortward of 10^ years. 

Key words: accretion — neutron star physics — Stars: binaries — X-rays:Stars 



1 INTRODUCTION 

X 1822—371 is a low mass X-ray binary and is the pro- 
totypical 'Accretion Disc Corona' (ADC) source (White & 
Holt 1982). The 5.57 hr orbital period of X 1822-371 (White 
et al. 1981) exhibits a quasi-sinusoidal variation in both the 
X-ray and optical (see also Mason & Cordova 1982; Hel- 
lier & Mason 1989). In addition, there is an approximately 
20 min long dip associated with partial obscuration of the X- 
ray source by the secondary mass-donating star (presumed 
to be filling its Roche lobe; see White et al. 1981, White & 
Holt 1982, Mason & Cordova 1982, Helher & Mason 1989, 
and Hellier et al. 1992). Thus X 1822-371 is believed to be 
a near edge on source; however, since the X-ray dip — hence- 
forth defined to be at zero orbital phase — is only partial, 
the X-rays are presumed to emanate from a very extended 
corona with a radius of the order of 3 x 10^" cm. 

Assuming a distance of 2 kpc (Mason & Cordova 1982a), 
the observed X-ray flux corresponds to an isotropic luminos- 
ity of Liao ~ 10^® ergs s"'^. The central X-ray source is ob- 



scured, however, and we only observe X-rays scattered into 
our line of sight. The intrinsic X-ray luminosity is therefore 
undoubtedly greater, perhaps substantially so (see §5). The 
properties and origin of this scattering corona are largely un- 
known. White & Holt (1982) postulated a corona, possibly 
optically thick, driven by a photoionizing radiation flux near 
the Eddington luminosity. In their model, a near Eddington 
luminosity is required to achieve the large scale height of 
the corona. Frank et al. (1987) postulated that the corona 
was due to interaction of the incoming accretion stream with 
the disc at the circularization radius. Viscous dissipation at 
radii within the circularization radius would in part lead to 
the vertical extent of the corona. 

The sinusoidal modulation is associated with obscura- 
tion by material with a vertical extent of order 10^*^ cm 
(White & Holt 1982; Hellier & Mason 1989). Frank et al. 
(1987) postulated that this obscuring material is also lo- 
cated at the disc circularization radius. However, the fact 
that the optical lightcurve also shows a sinusoidal modula- 
tion and a dip with approximately twice the duration of that 



© 0000 RAS 



2 S. Hemz & M. A. Nowak 



in the X-ray suggests that this obscuring rim is actually at 
twice this radius, i.e., closer to the disc tidal truncation ra- 
dius (Hellier & Mason 1989; see also §3 below). Again, this 
rim is associated with the interaction of the accretion stream 
with the accretion disc (Armitagc & Livio 1996; Armitage 
& Livio 1998). 

Perhaps the most confusing aspect of X 1822—371 has 
been its spectra. White et al. (1981) fit Einstein spectra with 
a flat, exponentially cutoff power law (photon index, V ~ 
— 1, cutoff energy ~ 17kcV). Furthermore, they required a 
broad (4 keV wide) Fc line component with equivalent width 
~ IkeV and a soft excess that they attributed to either a 
0.25 keV blackbody or possibly an Fe L complex. Hellier & 
Mason (1989), on the other hand, fit EXOSAT spectra with 
a fiat power law (P ~ —0.8), an ~ 2keV blackbody, and 
an Fe Ka line with 270 eV equivalent width. The flux of 
the blackbody had an implied emitting area consistent with 
1/400 of the surface area of a neutron star. Hellier & Mason 
(1989) therefore postulated that this emission is indeed from 
a neutron star surface and is scattered into our line of sight 
by a very optically thin corona. 

Several years later, Hellier et al. (1992) considered 
Ginga spectra of X 1822—371, and they attempted to fit 
the same model as for the EXOSAT data. Although such 
a model was the best simple fit that they could achieve, 
the fits were not adequate and yielded a reduced ~ 10- 
Furthermore, HcUior et al. noted that the dip at zero or- 
bital phase was fractionally larger at higher energies. They 
interpreted this fact as indicating that the corona was in 
fact optically thick. Hellier et al. (1990) claim that over the 
orbital timescales, the Fe Ka/f3 band shows greater variabil- 
ity than the continuum bands. More recently, Parmar et al. 
(2000) (hereafter POO) examined BeppoSAX and Advanced 
Satellite for Cosmology and Astrophysics (ASCA) data of 
X 1822—371. (We consider some of the same ASCA data in 
§2 below.) They consider a model comprised of Comptoniza- 
tion of a Wicn tail {kTw ~ 150 cV) in a very optically thick, 
low temperature corona (tcs ~ 23, fcTc ~ 7keV). Also re- 
quired in these fits are a blackbody with kT « 1.8 keV and 
an Fe Ka line with equivalent width 65-150 eV. This lat- 
ter feature possibly could be in reality two lines (Fe Ka//9), 
not adequately resolved from one another. POO found no ev- 
idence for either an Fe K-edge or O K-edge; however, they 
claim a detection of a 1.3-1.4 keV edge with t w 0.1-0.3 
that they associate with K-edges of Ne X and neutral Mg, 
or the L-edges of moderately ionized Fe. (See our discussion 
of §5, however.) 

The question of whether X 1822—371 contains a neutron 
star, black hole, or even a white dwarf primary also is yet 
unresolved. The lightcurves have been fitted with a model 
with a 1.4 Mq primary (Hellier & Mason 1989); however, 
as we discuss in §3 these fits do not uniquely determine 
the primary mass. To date, no variability has been detected 
that would uniquely point to a neutron star primary (see 
Hellier et al. 1990 for variability analysis of the Ginga data). 
Previous analyses only detected variability on the orbital 
time scales. 

Facing a new generation of X-ray satellites, we decided 
to revisit X 1822—371 with the help of predominantly un- 
published simultaneous archival Rossi X-ray Timing Ex- 
plorer (RXTE) and ASCA data. Our goal is to assess the 
evidence gathered so far and to point out what we believe 



Table 1. Observation Log. 
Satellite Instr. OBS ID Start Date Exp. Time 

RXTE PCA 10115-01 26 Sep. 96 21.6 kscc 
ASCA SIS,GIS 44015000 26 Sep. 96 28.5 ksec 
ASCA SIS,GIS 40019000 7 Oct. 93 38.5 ksec 

to be future avenues for cutting edge X-ray spectroscopy, as 
will be provided by the X-ray Multiple Mirror-Newton tele- 
scope (XMM-JVewton) and Chandra. First, we discuss the 
orbital ephemeris (§2.1) and the gross spectral variations 
over the orbit (§2.2). We then consider two separate spectral 
models which can be considered as broadly representing an 
'optically thick corona' and an 'optically thin corona' (§2.2). 
Using the RXTE lightcurves, we consider fits in multiple en- 
ergy bands (§3). Here we consider both optically thick and 
optically thin coronae, and furthermore we consider white 
dwarf, neutron star, and low mass black hole primaries. The 
high frequency (10"^-0.3 Hz) variability of the RXTE data 
is then considered (§4). We then discuss the implications of 
these analyses (§5) and present our conclusions (§6). 



2 DATA ANALYSIS 

In this paper we analyze a set of simultaneous RXTE and 
ASCA observations of X 1822-371 taken in 1996 and a 
separate set of ASCA observations from 1993, as summar 
rized in Table 1. Details of the data reduction procedure are 
summarized in the Appendix. For RXTE we only consider 
data from the Proportional Counter Array {PCA) , while for 
ASCA we consider data from both the Solid State Imaging 
Spectrometers {SIS) and the Gas Imaging Spectrometers 
{CIS). 

We use the simultaneous RXTE / ASCA data to per- 
form 'global', broad-band fits, and not for detailed line spec- 
troscopy. In what follows, we therefore have combined the 
two separate SIS spectra into one spectrum and likewise we 
have combined the two separate CIS spectra into a single 
spectrum. As we use the 1993 ASCA spectra for detailed 
line modelling, we do not combine the separate SIS spectra 
for those observations. 



2.1 Orbital Evolution 

Updated parameters on ephemeris and period changes have 
most recently been provided by POO, including all three data 
sets analyzed in this paper. To confirm the measurements by 
POO, we fitted the eclipse with a linear function attenuated 
by a Gaussian. Our fits agree with those of POO to within the 
error bars (indicated by parentheses below). POO et al. give 
a quadratic ephemeris for the eclipse midpoint as a function 
of orbital number, N, of 

Teci = 2445615.30964(15) + 0.232108785(50)N 

-1-2.06(23) X 10""N^ . (1) 

This implies a period derivative of 

rlP 

^ = 1.78(20) X 10"'° . (2) 
(See also Hellier et al. 1990.) Thus P~^P w 3.6 x 10® years. 
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Wc can now use this value to estimate the mass transfer 
rate in the system. We assume that a fraction 1 — / of the 
mass lost from the secondary is transferred to the primary, 
while the remaining fraction / is lost from the system at the 
location of the primary (e.g., in the form of a central disc 
wind or jot): Mi = — (1 — /)M2. The angular momentum loss 
is then j/J = /M2/M25^/(l + g) . If we define the usual 
mass ratio q = M2/M1, where Mi and M2 are the masses 
of the primary and the secondary respectively, we can write 
the mass accretion rate onto the primary as 



Ml = Ml 



3(1 - 9) P 



(1-/) 



1+i 



fq 



3 1 



For 

Ml 



( = 0.2 and Mi = 1.4 Mq (see §3) this gives 

= 2.1 X 10^^(1 -/)[l + 0.14/]-i gs-i 
= 3.3 X 10"* (1 - /) [1 + 0.14/]"^M(: " 



-1© yr 



(3) 



(4) 



Note that the Eddington accretion rate for a 1 Mq compact 
object is MEdd ~ 1-5 x 10^® (O.l/r?) gsec~^, where r? is the 
radiative efficiency. The second term in square brackets of 
eq. (3) is small for q < 0.7. If, on the other hand, most of 
the mass loss occurs in the form of a wind from the sec- 
ondary, the estimate for Mi is much less well-constrained. 
Since X 1822—371 is a LMXB, however, it is somewhat un- 
likely that the low-mass companion will have a very strong 
wind. 

The above mass loss rate implies a change of the Roche 
lobe radius of the secondary, R2, on comparably short time 
scales as the orbital period evolution. Specifically, one can 
show that given the above assumptions 



R2 



= -2 



5 fq 
(7 H — — 

6 ^^3(l + g) 



Ma 
M2 



(5) 



That is, for q ~ 0.2, the secondary's Roche lobe radius is 
expanding, even for large mass loss via a wind, and there- 
fore we expect strong mass transfer to be a short-lived phe- 
nomenon in this system. This conclusion is unaltered even 
if the (non-magnetic) wind mass loss occurs at radii as largo 
as the disc circularization radius. If instead we postulate 
conservative mass transfer but an unspecified source of an- 
gular momentum loss (perhaps from magnetic braking and 
gravitational radiation) sufficiently large to lead to R2 ~ 0, 
this implies an oven larger mass transfer rate than discussed 
above (Mi « 10^® g for a 1.4 M© primary), and also 
implies a characteristic angular momentum loss time scale 
of order 10^ years. We return to these considerations in §5. 



2.2 Spectral Analysis 

Given the uncertainties in the physical conditions in this 
system we felt it appropriate to keep spectral fitting to a 
phenomenological level. The fits we will present in the fol- 
lowing are meant to represent just the essential features of 
the two physical cases outlined in the introduction: optically 
thick and optically thin geometries. 

EXOSAT observations show an energy dependence of 
the modulation depth both for the sinusoidal component 
due to the accretion disc rim and the quasi-Gaussian eclipse 
due to the companion (HcUier & Mason 1989). Such a de- 
pendence might indicate a temperature stratification in the 
coronal region. We have tested the simultaneous RXTE and 
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Figure 1. Fractional depth of sinusoidal lightcurve variation (top 
panel, squares). The bottom panel (diamonds) show the product 
of amplitude times width of the Gaussian eclipse (both width and 
amplitude increase with energy). 



ASCA observations for such a trend by fitting a series of har- 
monics and a Gaussian to the X-ray lightcurves of the three 
instruments (PCA, SIS, and GIS). After inspection for con- 
vergence we truncated the fitting at the third harmonic. The 
fits confirm the trends reported by Hellier & Mason (1989). 
The fractional depth of the lowest order harmonic decreases 
with energy between 0.6 keV and 3 keV, and then rises again 
near the Fe Ka/K/3 region. The variation is small (of order 
10%), but it is statistically significant. Both the fractional 
depth of the Gaussian and its width increase with energy. 
The product of the two (proportional to the fraction of the 
emission blocked by the secondary) actually varies signifi- 
cantly with energy, roughly by 50%, which we plot in Fig.l. 

Due to the energy dependence of the modulation, some 
caution is in order when fitting global spectral models. We 
decided to divide the lightcurve into 5 parts (shown in Fig. 5) 
labelled A (the eclipse) through E. The detailed spectral fit- 
ting we have performed was carried out for each of these 
phases to test the dependence of different parameters on or- 
bital phase. We have also investigated the dependence of the 
PCA colours on phase. The strongest variations are found 
in the soft energy bands [2.5-3.7 keV/3.7-5 keV and 3.7-5 
keV/5-6.2 keV], where we find variations of order 30% to 
50%, with phase A being softest and phase D being hardest. 
In the more energetic bands [5-6.2 keV/6.2-10.2 keV and 
6.2-10.2 keV/10. 2-14.6 keV], the variation of colour with 
phase is weaker, i.e., of order 10%. These results are consis- 
tent with the analysis of BeppoSAX data by POO; however, 
the signal is not nearly as strong and as indicative as, for ex- 
ample, in the case of Her X-1. The strong phase-dependent 
behaviour in that case has been interpreted as being due 
to neutral hydrogen absorption and electron scattering in a 
large column depth medium (Stelzer et al. 1999). 

As a starting point we fitted test models to the com- 
plete data set (not divided into phases) . As with the previous 
analyses cited in the introduction, simple one-component 
models cannot produce a satisfactory fit — neither a single 
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blackbody, nor a thermal plasma, nor a cutoff power law 
adequately fit the data. Adding a second component can 
improve the fit significantly. We found that the combiner 
tion of a disc blackbody (Mitsuda et al. 1984) with either a 
simple blackbody or a cutoff power law (plus a Gaussian to 
model the iron Ka line) did produce marginally acceptable 
fits (residuals are plotted in Fig. 3). However, better fits can 
be achieved with physically more plausible models, which 
are consistent with the strong constraints that light curve 
fitting puts on the geometry of the system, which we will 
outline below. 

If the central object in X 1822—371 is in fact a neutron 
star, any thermal radiation from its surface is completely 
obscured by the accretion disc (HcUior & Mason 1989). The 
observed X-rays are produced either in the disc itself or in 
the corona. Lightcurve fitting indicates that the emission 
comes from a large extended source, and the (essentially re- 
quired) power law component in the spectra supports the 
notion that the corona is the main radaition source. We 
based our spectral fitting on this assumption and investi- 
gated two possible, though not exclusive, spectral models, 
one with an optically thick corona and one with an optically 
thin corona. 

The optically thick case: We expect the spectral signature of 
the hot corona to be of power law shape, with an exponential 
high energy cutoff produced by Comptonization, as is appro- 
priate for many LMXBs (White & Stella 1988). (Here a cut- 
off at low energies would indicate a relatively cool corona.) 
If the corona is optically thick, the spectrum will not contain 
any component originating from the central object, as such 
radiation will have been completely reprocessed. The accre- 
tion disk is nearly edge on and its atmosphere will partially 
cover the coronal emission. To account for this effect, we 
used a model consisting of a partially absorbed cutoff power 
law and an iron line. The partial absorption is modelled by 
co-adding an unabsorbed and an absorbed cutoff powerlaw 
with equal parameters. The normalization of the absorbed 
component, however, is multiplied by a constant, constcpi, 
which is related to the absorption fraction, /abs, by 



/a 



constcpi 



- const 



cpl 



(6) 



The absorption column in the partial absorption model is 
treated as a free parameter, and we include both electron 
scattering and neutral hydrogen absorption modelled us- 
ing the cross sections of Balucinska-Church & McCammon 
(1992). Our spectral model is, of course, a simplification be- 
cause in any natural situation a range of absorbing columns 
will occur. Higher levels of detail, however, are not warranted 
given our understanding of this source and the limits of the 
X-ray data. 

Having found a satisfactory overall spectral fit, we used 
the phase selected spectra to optimize the fits. Wc assumed 
that the physical parameters in the emitting region do not 
change from phase to phase and that the changes in the 
spectrum are entirely due to changes in obscuration of the 
corona both by the completely opaque disk and the partial 
absorption. We thus tied the powerlaw slope and the cutoff 
energy to be the same for all five phases, so that in the fitting 
process these parameters vary in unison. Similarly, we tied 
the iron line energy of all five phases together and, for lack of 
high spectral resolution from the PCA which dominates the 
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Figure 2. Phase C spectrum and best fit models in the optically 
thick case. Top panel: spectrum and model components. Dotted 
line: PCA components; dashed line: SIS components, (a) Unab- 
sorbed jjowerlaw: (b) absorbed powerlaw; (c) gaussian. From sec- 
ond panel down; residuals of different models. In sequential order 
downward; (1) disc blackbody plus cutoff power law plus narrow 
Gaussian at ~ 6.4 keV; (2) disc blackbody plus blackbody plus 
narrow Gaussian at ~ 6.4 keV; (3) partially absorbed cutoff pow- 
erlaw; (4) partially absorbed cutoff powerlaw plus iron line. The 
parameters for the fit in the bottom panel axe shown in Table 2. 



statistics in this region, fixed the iron line width to 0.1 keV. 
The foreground absorption column was also tied together 
for all five phases. Having combined the two SIS data sets 
and the two GIS data sets, we were left with a total of 15 
data sets, which we then fit simultaneously. The results of 
the joint fitting are presented in Table 2. The reduced is 

Xred = 0-78. 

It is worth noting that there is a correlation between the 
power law normalization (overall count rate), the internal 
absorption column, and the absorption fraction, with the 
internal absorption column being greatest in the brightest 
phases. Note also that the equivalent width of the iron line 
is very large (~ 250 eV, much larger than in the optically 
thin case). A more detailed interpretation of these results 
will follow in §5. 

The optically thin case: What distinguishes the optically 
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thin from the optically thick case is that radiation from the 
central source can be scattered into the observer's line of 
sight by the corona, without having its input spectrum be- 
ing completely distorted. One might therefore expect to see 
a blackbody component in combination with a cutoff power 
law. In principle, both of these components will be partially 
absorbed, as Eissumed in the optically thick case. Since the 
opically thick model already fits the data satisfactorly, how- 
ever, such a model would be an over-parametrization given 
the limits of the data. To highlight the differences between 
the two models, we therefore neglected partial absorption in 
our fits that included blackbody emission. 

Such a model can fit the data reasonably well, leaving 
residuals at both the low and high energy ends. The low 
energy residuals could be caused by a line complex, as also 
would be expected to emanate from an optically thin (prob- 
ably photoionized) corona. We have modelled these low en- 
ergy residuals with the meka plasma model; (Kaastra 1992, 
and references therein). Once again, we fitted the five phase 
selected spectra simultaneously, and we constrained the tem- 
peratures of the blackbody and the thermal plasma to be 
uniform throughout all five phases. We also constrained the 
line energy and both the power law index and cutoff energy 
to be uniform throughout the phases. The resulting best fit 
with x^red = 0-^4 is presented in Table 3 and in Fig. 3. 

We note that the lack of absorption in the best fit 
model is somewhat inconsistent with the predicted GOLDEN 
NH value. Forcing the hydrogen column to the GOLDEN value 
worsens the fit slightly, loading to residuals below 0.8 keV. 
Since the exact distance to X 1822—371 is unknown and the 
lack of partial absorption in the optically thin case is an 
idealization anyway, we decided to leave the clarification of 
this issue to Chandra and XMM observations. 

In both cases, the high energy end of the spectrum 
shows some systematic residuals, present also in the fits 
based on disc emission models. We suspect that at least 
some of this effect can be attributed to systematic errors in 
the RXTE response, in particular, differences in the high en- 
ergy slopes of ASCA and RXTE. However, the Ginga data 
reported by Hellier et al. (1992) also seem to show this hard 
excess, which would argue against a purely instrumental ef- 
fect. In this paper we will not attempt a physical explanation 
of this feature. The failure of Astro-E with its high energy 
capabilities is particularly unfortunate in this respect. 

We then tried to determine the structure of the iron 
line. Since the sensitivity of the SIS has been deteriorating, 
we decided to use only the 1993 SIS datasct from the 1993 
ASCA observation. We limited the fits to the energy range 
from 3-10 keV and approximated the continuum spectrum 
around the line by a cutoff power law. Two narrow lines pro- 
duce the most significant improvement in , from 587/458 
to 421/454, in agreement with POO. The line energies are 
6.4 keV and 7keV, consistent with a cold Kq line and either 
K/3 or ionized Ka. The equivalent widths of the two lines 
are on the order of 80 eV and 40 eV respectively. 



3 MODELLING THE LIGHTCURVE 

Models of the X 1822-371 lightcurve, for both the opti- 
cal and X-ray energy bands, have been presented by White 
& Holt (1982), Mason & Cordova (1982a), and Hellier & 



Model: const. ■ phabs(NHphabs) [bbody(Tbbd, Nbbd) + cutoffpl(rcpi, Ecutcpi, Ncpi) + gauss(Egss, o-gss, Ngss) + raymond(Tray , 1, 0, Nray)] X^/dof = 



3452 
4665 



: 0.74 



Table 3. Best fit parameters for the 'optically thin model'. Physical parameters (blackbody temperature, power law slope and cutoff, plasma temperature, and iron line energy) 
are tied to stay constant as a function of phase, denoted by the symbol " . The reduced for this fit is x^red ~ 0-74, with 4655 degrees of freedom. We have frozen the iron 
line width to 0.1 keV. Similarly, redshift and metallicity of the Raymond-Smith component were frozen to their canonical values (0 and 1 respectively). 
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disc rim profile is then linearly extrapolated at intermediate 
phases. Fit parameters for these models have included the 
radius of the X-ray emitting corona, usually assumed spher- 
ical for simplicity, the heights and phases of the rim nodes, 
the mass of the compact object, and the inclination angle 
of the system with respect to our line of sight (see Hellier & 
Mason 1989). 

Here we elaborate upon these models in two ways. First, 
due to the excellent statistics of RXTE, we are able to di- 
vide the RXTE lightcurves into five distinct energy bands 
covering PC A pha channels 8-11, 12-15, 16-19, 20-27, and 
28-32 (2.9-4.4 keV, 4.4-5.8 keV, 5.8-7.2 keV, 7.2-10.2 keV, 
and 10.2-12.0 keV, respectively). We fit each of these energy 
band lightcurves simultaneously. Second, we model both an 
optically thin and an optically thick spherical corona. Specif- 
ically, we model the X-ray emission via a uniform emissivity 
or a uniform surface brightness. We do not expect (and as 
we confirm below; see also Hellier & Mason 1989) strong 
differences between the two cases. Given a completely edge- 
on viewing angle, for a given disc rim height the fraction of 
emission that is obscured and the emission-weighted mean 
height of the unobscured radiation are nearly identical for 
both uniform surface brightness and uniform emissivity. In 
our fits, we allow both the overall fiux normalization and 
the radius of the sphere to vary for each energy band; how- 
ever, the disc rim parameters (radius of the disc, -Rrim, node 
heights, H , and node phases, 0) are held fixed. Additional 
fit parameters are the mass of the compact object (Mi) and 
the inclination of the system with respect to our line of sight 
(i). 

In part these parameters are determined by utilizing 
the fact that X 1822—371 has a known binary period of 
P = 5.57 hr and a measured projected primary orbital ve- 
locity of Ki = 70 kms~^ (Mason et al. 1982; Cowley et al. 
1982). We highlight the systematic uncertainties in our fit 
parameters by varying this latter parameter by -|-30 kms~^ 
(see Cowley et al. 1982). Furthermore, we search for model 
fits with Ml ~ 0.3, 1.4 and 2.5 Mq. That is, we fit mod- 
els for a 'white dwarf, 'canonical neutron star', and 'low 
mass black hole' primary. The latter primary mass yields a 
secondary mass consistent with Roche lobe overflow from a 
main sequence star (see below). 

In order to better interpret the results of our fits, we 
first present estimates of the characteristic accretion sys- 
tem size and mass scales. Defining the secondary to primary 
mass ratio, q = M2/M1, Newton's laws give q ~ 0.21 for 
Ml = 1.4 Mq, Ki = 70 kms"\ and i = 85°. Henceforth, 
these parameter values shall be referred to as the 'nominal 
parameters'. As noted by Mason et al. (1982), this mass 
ratio is roughly half that expected from the mass-period re- 
lationship for a lower main sequence secondary. If instead we 
choose Ml = 2.5 Mq, ATi = 100 km s"\ and i = 85°, one 
obtains a secondary mass of M2 = 0.62 Mq , as one would ex- 
pect for the period-mass relationship of a main sequence star 
(Prank et al. 1992). The above system parameters yield bi- 
nary separations of a ~ 1.3-1.6 x 10^^ cm, respectively. Over 
the parameter ranges of interest to us, q varies roughly as 
(/fi/sini)Mj and the binary separation weakly varies 
as (l + g)i/^ M'/\ 

The disc circularization radius can be determined via 
the approximation 



R^i,^ = {1 + q){0.5 - 0.227 log g)^ a (7) 

(Frank et al. 1992), where for the nominal parameters 
Rcirc ~ 0.22a « 2.9 x 10^" cm. As a fraction of the binary 
separation, this ratio is fairly constant for the parameters of 
concern to us. The disc tidal truncation radius is approxi- 
mately given by Rt « 0.9 Ri, where Ri is the Roche lobe 
radius of the primary. Approximating this radius (Eggleton 
1983) as 

Ri ^ 0-49g-"/^ 
a ~ 0.6g-2/3-^ln(l-hg-i/3) 

yields Rt « 0.47 o 6.1 x 10^^° cm for the nominal parame- 
ters. Again, expressed as a fraction of the binary separation, 
this radius is relatively constant for the parameters of inter- 
est to us. 

Equation 8, with q replaced by also yields Tig /a 
(~ 0.22 for the nominal parameters), where R2 is the Roche 
lobe radius of the secondary. For a near edge on system 
with a brief (relative to the orbital period) eclipse, 7?2/a 
is approximately proportional to the fraction of the binary 
period that can be eclipsed by the secondary. This fraction 
is roughly proportional to q^^^ oc M^^^^{Ki/ sini)^^^. As 
the duration of the eclipse is a fixed fraction of the binary 
orbital period, we therefore expect the coronal radius, ex- 
pressed as a fraction of the binary separation, to be weakly 
dependent upon Ki / sin i and almost completely indepen- 
dent of Ml (see also Mason & Cordova 1982). As we show 
below, we indeed can fit a relatively wide range for Mi*. 
The mass ratio, q, however, is more tightly constrained by 
the lightcurve fits. 

We divide the lightcurves into 50 phase bins in each en- 
ergy band. In lieu of statistical errors, which were very small 
given the large effective area of RXTE, a systematic error 
of 4.2% was added to each phase bin. This represented the 
average variance of the lightcurve from orbit to orbit over 
the X 1822—371 period. (Approximately four orbital peri- 
ods were measured in part or in whole.) For the optically 
thin case, the corona was divided into 6 evenly spaced ra- 
dial zones, 50 evenly spaced zones in the azimuthal angle, 
(j), and 50 zones evenly spaced in ij, = cos 9 in the polax an- 
gular direction, 0. The optically thick corona had 50 x 50 
evenly spaced zones in 4> and p. We assumed a disc rim 
with 9x2 nodes symmetrically placed above and below 
the disc mid-plane. Ray tracing was performed to deter- 
mine whether emission from a given coronal element inter- 
cepted the disc rim. Near zero phase, the blink subroutine 
from Keith Home's eclipse mapping code (Horne 1985) was 
used to determine if the secondary was blocking the line of 
sight to the emission element. In all, there were 31 fit pa- 
rameters (primary mass, system inclination, 5 coronal radii 
and emissivities/surface brightnesses, disc rim radius, and 9 
node heights and phases). For the fits, all energy bands were 
weighted equally, and the amoeba subroutine from Press et 
al. (1992) was used in the minimization. Fit results are 
presented in Tables 4 and 5. 

Fig. 5 shows the best fit optically thick corona 
model for the 4.4-5.8 keV energy band, as well as 

* If we insist that the (main sequence) radius of the secondary 
not be greater than R2, then we require Mi ^ 2.5. Hence we only 
consider primary masses of fs 0.3, 1.4, and 2.5 Mq. 
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Table 4. System and Coronal I'arainctcrs for Fits to the X 1822—371 RXTE Lightcurves. All length scales are in units of the binary 
separation, a. Rc^j, is the radius of the optically thin or thick, uniformly emitting corona for RXTE pha channels a-b. Rrim is the radius 
of the obscuring rim. Mi is the mass of the central compact object, and i is the angle of the normal to the binary orbital plane with 
respect to our lino of si,5^ht. Values ,i^i\'en are assuiiiiiiji a projeeted primary veloeit\' of A'[. 



Model 


Ml 
(Mo) 


(km s-i) 


i 

n 




^C12_15 




■^C20_27 


■^C28_32 




xVdof 


Thin 


0.32 


70 


80.7 


0.220 


0.216 


0.210 


0.211 


0.211 


0.38 


188/219 


Thin 


1.40 


70 


83.4 


0.219 


0.217 


0.216 


0.216 


0.215 


0.41 


154/219 


Thin 


2.47 


100 


82.5 


0.220 


0.216 


0.213 


0.213 


0.212 


0.40 


175/219 


Thick 


0.31 


70 


81.1 


0.216 


0.213 


0.209 


0.210 


0.210 


0.40 


181/219 


Thick 


1.38 


70 


83.7 


0.215 


0.214 


0.213 


0.213 


0.213 


0.39 


180/219 


Thick 


2.48 


100 


83.0 


0.216 


0.213 


0.212 


0.212 


0.212 


0.41 


180/219 



Table 5. Disc Rim Parameters for Fits to the X 1822—371 RXTE Lightcurves. All length scales are in units of the binary separation, 
a, and phases are in units of the binary phase, rp and H are, respectively, the phase and height of disc rim nodes Ni-Ng. Rim heights 
inbetween the nodes are determined by linear interpolation. Values given are assuming a projected primary velocity of Ki. 

Model Ml Ki Ni N2 N3 N4 N5 Ne N7 Na Ng 

(Mo) (kms-i) 



Thin 


0.32 


70 


4> = 


0.111 


0.295 


0.321 


0.408 


0.526 


0.750 


0.813 


0.934 


0.943 








H = 


0.110 


0.093 


0.075 


0.102 


0.099 


0.136 


0.132 


0.109 


0.108 


Thin 


1.40 


70 


</> = 


0.107 


0.202 


0.317 


0.399 


0.564 


0.741 


0.846 


0.882 


0.938 








H = 


0.104 


0.105 


0.080 


0.097 


0.103 


0.129 


0.132 


0.096 


0.117 


Thin 


2.47 


100 


4> = 


0.124 


0.299 


0.317 


0.413 


0.537 


0.768 


0.830 


0.889 


0.907 








H = 


0.106 


0.090 


0.078 


0.099 


0.100 


0.133 


0.132 


0.100 


0.112 



Thick 


0.31 


70 


4,= 


0.112 


0.291 


0.323 


0.410 


0.535 


0.750 


0.831 


0.924 


0.928 








H = 


0.111 


0.090 


0.074 


0.101 


0.100 


0.139 


0.134 


0.103 


0.108 


Thick 


1.38 


70 


</,= 


0.103 


0.292 


0.359 


0.441 


0.553 


0.756 


0.828 


0.896 


0.911 








H = 


0.111 


0.088 


0.082 


0.103 


0.101 


0.137 


0.132 


0.104 


0.112 


Thick 


2.48 


100 


4 = 


0.111 


0.296 


0.343 


0.419 


0.537 


0.765 


0.819 


0.922 


0.930 








H = 


0.109 


0.088 


0.079 


0.100 


0.101 


0.137 


0.132 


0.104 


0.110 



the fitted disc rim profiles for both the optically 
thick and thin models, both for Mi ~ 1.4 Mq, 
Ki = 70 km s~^. (Animations of the fits can be viewed 
at http: //rocinante . Colorado . edu/~ heinzs/1822/ . ) Wc 
note that using fewer coronal grid points did not adequately 
resolve the coronal emission; however, numerical discrete- 
ness in the fitting process led to uncertainties of Ax^ ~ 2 
in any fits with a greater number of coronal grid points. In 
Table 4 we therefore only present the systematic errors asso- 
ciated with varying the projected primary velocity, Ki, be- 
tween 70 and 100 km s~^ and from searching for minima 
near primary masses of 0.3, 1.4, and 2.5 Mq. We also experi- 
mented with varying the number of nodes in the fit; however, 
9 nodes were sufficient to produce a reduced ~ 1- 

Both the optically thick and optically thin coronal mod- 
els fit the data nearly equally well, with ranging from 
154-188 for 219 degrees of freedom. Although there is a 
slight preference for an optically thin corona with Mi « 
1.4 Mo, given the simplicity of our assumptions we do not 
consider these differences to be strongly significant. 

Most of the fits showed a weak trend for the coronal 
radius to decrease with energy, as would be expected from 



the energy-dependent fractional modulation shown in Fig. 1. 
The fractional change in the coronal radii was as large as 5%; 
however, this is still somewhat smaller than the ~ 10% found 

in the two energy channel (l-6keV and 6-30 keV) Ginga 
data fits of Hellier et al. (1992). The other size scale param- 
eters are consistent with previous models of the lightcurve 
(White & Holt 1982; Mason & Cordova 1982a; Hellier & Ma- 
son 1989). Specifically, the disc rim height shows a inininmm 
near phase 0.3-0.35, and a maximum near phase 0.75-0.85. 
The disc rim height is fs 0.1a w 1.5 x 10^^° cm. The radius of 
the disc rim is consistent with being at the tidal truncation 
radius, which is reasonable if the disc rim is due to the in- 
teraction of the secondary's accretion stream with the disc 
(see, for example, Armitage & Livio 1996, 1998, and refer- 
ences therein). The coronal radii are all « 0.2o « 3x10^^° cm, 
i.e. consistent with the corona extending all the way to the 
disc circularization radius. 

We find a system inclination of i w 81°-84°, which 
is also consistent with previous models of the lightcurve 
(Hellier & Mason 1989). As expected from the analysis dis- 
cussed above, the fits for Mi « 0.3, 1.4, and 2.5 Mo were 
nearly identical. Thus, based upon the lightcurve analysis 
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Figure 5. Top panel: Fblded X 1822-371 lightcurve in RXTE 
pha channels 12-15 (4.4-5.8 keV) plus the best fit model (solid 
line) for a uniformly emitting, optically thick corona obscured by 
a disc rim. (Dotted lines and labels A-E indicate the regions of 
the phase-resolved spectral fits.) Middle panel: The best fit disc 
rim profile, in units of the binary separation. The dashed lines 
show the projected positions (relative to the disc rim, with re- 
spect to our line of sight) of the coronal poles. Bottom panel: The 
same as the middle panel, but instead for the best fit optically thin 
corona model. Note: Both of the above models assume a projected 
primary velocity of K\ = 70 km and Mi ss 1.4 Mq. Anima- 
tions of the X 1822—371 system and lightcurve can be found at 
http : / /rocinante . Colorado . edu/~ heinzs/1822/ . 



alone, it is impossible to determine whether X 1822—371 
is a white dwarf, neutron star, or a black hole. Further 

observations, especially any that can independently mea- 
sure the velocity of the secondary, are required to break 
the degeneracy in the primary mass values. We note that 
to this end, Haxlaftis et al. (1997) place a lower limit on 
the companion velocity of K2 ^ 225 km s~^. If the lower 
end of this value is accurate and Ki — 70 km s^^, then 
q = K1/K2 ~ 0.3 (consistent with the above estimates) and 
Ml Ma « 0.6 Mq oc (/sTi + , i.e., X 1822-371 is a 



white dwaxf system. 



4 VARIABILITY ANALYSIS 

Extremely high time resolution data were available from 
the RXTE observations. Accordingly, we created lightcurves 

with 2^^''s resolution in order to search for high-frequency 
variability. No significant variability was detected at fre- 



with 0.5 s resolution. We chose three energy bands cover- 
ing PCA pha channels 5-10, 16-19, and 28-41 (1.8-4 keV, 
5.8-7.2 keV, and 10.2-15.3 keV respectively). These energy 
bands were chosen to represent a low energy band, an Fe 
Ka/P band, and a high energy band. Energies higher than 
15keV were background dominated. 
We consider three measures of the variability: the power 
spectral density (PSD) in each energy band, and the Fourier 
froquoncy-dcpondont time lags and coherence function be- 
tween variability in different energy bands. A discussion of 
Fourier techniques in specific, and timing analysis in gen- 
eral, has been presented by van der Klis (1989). Here we 
apply these Fourier analysis techniques in the same manner 
as for our RXTE observations of Cyg X 1 (Nowak ct al. 
1999). Specifically, we used the same techniques for esti- 
mating: deadtime corrections (Zhang et al. 1995; Zhang & 
Jahoda 1996); the error bars and Poisson noise levels of the 
PSD (Leahy et al. 1983; van der Khs 1989); the error bars 
and noise levels for the coherence function (Bcndat & Pier- 
sol 1986; Vaughan & Nowak 1997); and the error bars and 
noise levels for the Fourier frequency-dependent time lag 
between hard and soft photon variability (Bendat & Piersol 
1986; Nowak et al. 1999). Results of these analyses for the 
lowest and highest energy bands are presented in Fig. 6. In 
this figure, the PSDs are normalized according to Belloni & 
Hasinger (1990) wherein integrating over Fourier frequency 
yields the moan square variability relative to the square of 
the mean of the lightcurve. 

All three PSDs had roughly comparable shapes and am- 
plitudes. Specifically, they show evidence for a rise in power 
at frequencies ^ 3 x 10~^ Hz, broad peaks at / ~ 0.02 and 
0.08 Hz, and reasonably sharp declines for / ^ 0.1 Hz. No 
power in excess of noise is seen above 0.3 Hz. The root 
mean square (rms) variability increases from soft to hard 
energies. From softest to hardest energy bands, the rms vari- 
abilities are 8%/7%, 10%/8%, and ll%/9% over the lO'^^- 
0.3Hz/10^'^ 0.3 Hz range. The up-turn at low frequency is 
most likely due to variability associated with the orbital time 
scales. (Due to low count rates, phase resolved variability 
studies are very difficult.) The coherence at / ^ 10~^ Hz be- 
ing slightly lower than that at higher frequencies might be 
related to an admixture of intrinsic variability of the corona 
with variability associated with the orbital period, although 
overall the varibiality between the highest and lowest energy 
bands is well-correlated. We note that 7 9% rms variabil- 
ity in the 10^^-0.3 Hz range is consistent with the variabil- 
ity seen in both low (< 5% I/Edd) and high (> 30% Lfidd) 
fractional Eddington luminosity neutron star and black hole 
sources. 

The lack of high frequency variability may be due to 
scattering over large distances and/or optical depths (see 
Nowak & Vaughan 1996, and references therein). Specifi- 
cally, if the variable lightcurve is first passed through a scat- 
tering medium of optical depth t and size D, one expects a 
cutoff in the power spectrum at a frequency, /cut, given by 



27r/c, 



c c 



(9) 



quencies ^ 1 Hz; therefore, we created a series of lightcurves Thus, given Z) ~ J?c ~ 0.2o ~ 3 x 10^° cm, the cutoff fre- 
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quency should be /cut ~ 0.2 min[T~'^, 1] Hz^, consistent with 
what is observed here. 

Scattering may also lead to the observed variability lags. 

Over the narrow frequency region for which one can obtain 
lag measurements, the softest energy band is seen to lag the 
hardest energy band by w 1 s. (Note, the statistics wore not 
sufficient to measure accurately the time lag of the variabil- 
ity in the middle energy band with respect to that in either 
of the other two energy bands.) Such a lag is consistent with 
a central X-ray source being reprocessed downward in en- 
ergy by scattering within a corona of radius/scattering path 
length of Oder 3 x 10^^° cm, consistent with the fits of §3. 



5 DISCUSSION 

Having laid out the observational facts, we now turn to their 
interpretation. Several key issues about X 1822—371 have 
surfaced in the earlier sections, which we will now address, 
along with some open questions which relate to earlier stud- 
ies of this object. 



5.1 Detailed Spectral Fitting 

As we have shown, several fundamentally different models 
can explain the ASCA and RXTE spectra almost equally 
well. While formally satisfactory, the spectra we presented 
in §2.2 still show some systematic residues, which could be 
removed by addition of features like absorption edges and 
lines. For example, POO argue that there is clear evidence 
for an edge around 1.3 keV. We do see a hint of such an ef- 
fect, especially in the optically thick case, but it is not clear 
if the origin is really an absorption edge. The same feature 
might be produced by an emission line complex at 1 keV 
(as demonstrated in the optically thin case, where the ad- 
dition of a Raymond-Smith plasma removed the systematic 
residues at the low energy end). We also note that the edge 
in the fits by POO falls very close to the position where two 
major spectral components (blackbody and Comptonization 
spectrum) cross in their fits. 

Similarly, the structure of the iron line region is all but 
clear. As remarked by POO, the ratio of the two narrow lines 
we fit to the data is off firom the theoretical value expected 
for Ka to K/3, but ionized line emission could explain this 
fact. In the optically thin models the ionization parameter in 
the corona is much larger than one, which would strengthen 
such an argument. Both of these questions will undoubt- 
edly be answered by upcoming XMM-Newton and Chandra 
observations. 



t For an input spectrum confined to a narrow energy band scat- 
tered into a narrowly observed output energy band, the cutoff 
frequency can be given by 27r/cut ~ c/D, even for r > 1. This 
is because for such a situation, fixing the energy band of tfic ob- 
served output is essentially fixing the number of scatters that 
the observed photons have undergone, and therefore narrows the 
dispersion in the photon arrival times. See Nowak & Vaughan 
(1996), and references therein. 



5.2 Optically Thick vs. Thin 

Both in spectral modelling and in fitting the light curve, we 
used two general models, based on whether the corona is 
optically thick or thin. We found that both light curve and 
spectrum are degenerate with respect to this distinction, and 
optically thick and thin models reproduce these data equally 
well. We will now discuss the physical interpretation of the 
two cases (see also §5.4). 

• In the optically thick case we modelled the spectrum 
as a partially absorbed power law, produced by reprocessed 
radiation from the corona, passing through a column of cold 
gas (the atmosphere of the outer disc). From the spectral fit 
we can then deduce the covering fraction and optical depth 
of the partial absorber. These numbers are given in Table 2. 

We can understand the above mentioned correlation be- 
tween the partial absorption column, the absorption frac- 
tion, and the overall fiux through the following picture: the 
X-rays originate from an approximately spherical corona (at 
uniform surface brightness). The accretion disc rim is com- 
pletely opaque and responsible for most of the modulation 
in the light curve, due to differences in its height (and thus 
its covering fraction). Above the accretion disc, an atmo- 
sphere of cold gas absorbs/scatters a fraction of the light 
from the corona behind it. The geometric covering fraction 
would roughly be given by the ratio of the atmosphere's 
scale height to the coronal radius. As the rim height changes 
through the binary phases, the covering fraction of the par- 
tial absorber also changes (see Table 2). 

Higher latitudes of the corona are relatively less absorbed, 
thus for binary phases where the the disc rim is higher, the 
ratio of unabsorbed to absorbed coronal flux will be greater. 
At the same time, the average column depth of the absorber 
should be smaller, since the gas density is lower at higher 
latitudes. The trend in Fig. 1 can also be understood this 
way: the absorbed component of the power law contributes 
primarily hard X-rays to the spectrum. Since this part of 
the spectrum originates at low latitudes, it is most affected 
by obscuration, thus the larger modulation depths at high 
energies. The implied scale heights and column depths of 
Table 2 are consistent with what would be expected from an 
X-ray heated atmosphere above the disc between the edge 
of the corona and the disc rim (White & Holt 1982). 

The fact that we only see a small fraction of the luminos- 
ity inferred from the orbital evolution can be explained if the 
corona is optically thick to scattering while relatively opti- 
cally thin to absorption. In that case it will act as a mirror 
and transmit only a fraction {1 + t)~^ of the incident radia- 
tion. For a ratio of observed to inferred luminosity of ~ 200, 
modulo a geometric covering factor, this would imply an 
optical depth of r 200 (which justifies the assumption of 
an optically thick corona in this model). As we will discuss 
below (see §5.4), a large fraction of the energy could also 
leave the system in a wind (Adiabatic Inflow-Outflow Solu- 
tion, ADIOS; Blandford & Begelman 1999) or, in the case 
of a black hole primary, disappear down the horizon (op- 
tically thick Advection Dominated Accretion Flow, ADAF; 
Narayan & Yi 1994). 

Table 2 shows that the equivalent widths of the iron lines 
are rather large, of order 250 eV. A simple estimate of the 
expected EW produced solely within the partial absorber 
by the incident coronal spectrum is of the order of 50 eV. 
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Figure 6. Left: Power spectral density times Fourier frequency in the 1.8-4 keV and 10.2-15.3 keV energy bands (filled and clear 
diamonds, respectively). The solid line (1.8-^ keV) and dashed line (10.2-15.3 keV) show the residual noise levels, which correspond to 
the expected amplitudes of positive 1-<t fluctuations above the mean value of the Poisson noise PSDs. Centre: Coherence function between 
the variability in the 1.8-4 keV and 10.2-15.3keV energy bands. Right: Fourier frequency-dependent time lags between variability in the 
1.8-4 keV and 10.2-15.3 keV energy bands. Filled diamonds indicate the soft band lagging the hard band, and clear diamonds indicate 
the hard band lagging the soft band. The solid line corresponds to the magnitude of the expected 1-<t Poisson noise level. 



Uncertainties in the metallicity or the geometry of the sys- 
tem (i.e., if the partial absorber is subject to a larger flux 
than seen by the observer) might account for such a discrep- 
ancy. Note also that the partial covering model introduces 
an iron absorption edge at 7.1 keV, which is leading in part 
to the large fitted equivalent widths. As noted by POO and 
references therein, the line ratios are still problematic if the 
higher energy line is an Fe K/3 line, as opposed to ionized 
Fe Ka. Models of resonant scattering of the Fe Ka line by 
the partial absorber in the outer disc rim region may help 
explain some of this discrepancy. Finally, we note that the 
optically thick corona itself could well be the origin of strong 
iron line emission. 

• The key question for the optically thin case is whether 
or not there is unambiguous evidence for blackbody emis- 
sion from the surface of the compact object that is scat- 
tered into our line of sight. Related to this is the question 
of whether or not there is clear evidence for coronal line 
emission, i.e. a thermal plasma component. As the optically 
thick model, with partial absorption, fits the spectrum ade- 
quately, evidence for these components cannot be claimed to 
be unambiguous. Ignoring the partial covering (due to the 
limits of our data and models), however, these components 
are seen to be allowed in the fits, thus there absence cannot 
be definitively shown either. 

Assuming the optically thin scenario to be correct, the X- 
ray emission is again produced in the corona. The blackbody 
part of the continuum is produced cither on the neutron star 
surface (if the central object is in fact a neutron star) or the 
inner disc, and is then scattered into our line of sight by 
the corona. The power law part of the spectrum stems from 
radiation that is Comptonized in the corona. If the corona 
is optically thin, the fraction of the total X-ray luminosity 
that reaches the observer is « t. The optical depth of the 
corona would thus be r ^ 1/200, once again modulo the 
geometric covering factor. We note again that a fraction of 
the accretion energy might be carried away in a wind or 
advectcd (in the case of a black hole primary). 

Note that the modulation seen in the power law is stronger 
than the modulation seen in the blackbody. We can interpret 



this dependence geometrically if we further postulate that 
the power law emission is produced at lower coronal lati- 
tudes, in which case a larger fraction of it would be subject 
to obscuration both by the rim and the companion. The vari- 
ations shown in Fig. 1 would be natural in such a scenario, 
since the black body component emanates from the whole 
corona, is relatively less modulated, and contributes most 
strongly at low and intermediate energies. Again, however, 
a partial covering model may be applicable for the optically 
thin case as well as for the optically thick case. 

The fitted iron line in this case is weaker than for the opti- 
cally thick case, but due to the neglect of partial absorption 
this model does not contain an intrinsic absorption feature 
around 7 keV. AH other remarks with respect to the line 
strength in the optically thick case hold here too. 

Both these models can explain the basic features of light 
curve and spectrum, we therefore have no way to chose one 
over the other. In fact, reality might fall in between these 
two cases. One might imagine a model where the optical 
depth of the corona decreases with height. We would thus 
see scattered blackbody radiation from high latitudes and a 
partially absorbed power law from the optically thick parts 
closer to the disc. Hopefully, future X-ray and optical obser- 
vations will resolve this degeneracy. 



5.3 NS, BH, or WD? 

Lacking detections of X-ray pulses or Type I X-ray bursts 
that would unambiguously point toward a neutron star pri- 
mary, we must rely on more indirect arguments for deter- 
mining its nature. As discussed in §3, the fits to the X-ray 
lightcurves arc ambiguous in this regard. Of the possibil- 
ities, we find the white dwarf primary scenario the least 
likely for two principle reasons. First is the fact that we ob- 
serve am X-ray flux of order 10^® ergs a~^. Eq. 4, assuming no 
mass loss and a more typical white dwarf radiative efficiency 
of Tj ~ 10^'', would lead to a total accretion luminosity of 
« 4 X 10^^ ergs s~^. Although by assuming a white dwarf 
primary we decrease the inferred size of the disc and thereby 
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decrease the inferred distance (Mason & Cordova 1982a) and 
luminosity by a factor of ~ 0.6 and ~ 0.4, respectively, we 
would still require that we are viewing of order 100% of 
the accretion luminosity despite the presence of a scattering 
corona and the near edge-on inclination of the system. Sec- 
ond is the implied mass of the secondary, M2 ~ 0.06 Mq, 
which is slightly small given the then inferred Roche lobe 
radius of the secondary, J?2 « 2 x 10^° cm fa 0.3 R©. Again 
we note, however, that if the lower values of Ki = 70 km s~^ 
and K2 — 225 km s^^ arc correct (Harlaftis ct al. 1997), than 
a white dwarf primary is the preferred model. 

Of the oter two possibilities, neutron star or black hole, 
there is little to distinguish between them. Given the edge-on 
inclination and the implied large scattering ijaths, it is not 
surprising that we do not detect any variability associated 
with a neutron star spin period, as discussed in §4. Further- 
more, given little or no mass loss, the implied accretion rate 
onto the neutron star would be suflScient to suppress any 
Type 1 bursting behaviour (Bildstcn 1995, and references 
therein). We therefore expect a neutron star and low mass 
black hole to look nearly identical. Distinguishing between 
these two scenarios, therefore, will require a combination of 
more careful observations of the secondary and identifying 
a plausible evolutionary scenario for this system. 

5.4 Winds, Advection Domination, and the 
Evolutionary History of X 1822-371 

As discussed in §2.1, the implied mass transfer time scale 
is Tmt ~ lO^years, which is compatible with mass transfer 
on a thermal time scale (see Kalogera & Webbink 1996, and 
references therein). One interesting possibility is that we are 
viewing the X 1822—371 system towards the end of a ther- 
mal time scale mass transfer phase (Kalogera & Webbink 
1996). Specific scenarios have been discussed for Cyg X-2 
(King & Ritter 1999; King & Begelman 1999) and SS 433 
(King et al. 2000). At the onset of mass transfer the system 
begins with a mass ratio g > 1. As the secondary loses mass, 
both the binary separation and the secondary Roche lobe 
shrink. A common envelope phase is avoided by strong mass 
loss from the primary, perhaps in an advection dominated 
accretion phase cliaracterizcd by a strong wind emanating 
from large radii (King & Begelman 1999; King et al. 2000). 
Even as q decreases below unity, mass transfer can continue 
to be driven on the thermal time scale. For the Cyg X-2 
system (known to be a neutron star primary; Smale 1998), 
it is hypothesised that the system began with an « 3.5 Mq 
secondary, but now shows an « 0.5 Mq secondary with a 
radius w 7 R© (King & Ritter 1999). 

The above mentioned systems, however, cannot be ex- 
actly analogous to that of the X 1822—371 system. Cyg X- 
2, for example, has a 9.84 day orbital period and a sec- 
ondary luminosity of « 150 Lq, both far larger than in 
X 1822—371. In the optical bands, aside from X-ray heat- 
ing of the secondary, it has traditionally been assumed that 
the X 1822—371 emission is dominated by the accretion 
disc and disc rim, with the total optical luminosity being 
iopt ~ 10 L0 (Mason & Cordova 1982a). Likewise, models 
of the UV emission have been presumed to be dominated 
by the disc and the corona (Mason & Cordova 1982b), al- 
though perhaps the far UV allows the greatest room for 
significant contributions from the secondary. In addition. 



obtaining the currently observed 5.57 hr orbital period via 
mass transfer alone requires that during the earlier epoch 
wherein 0.2 ^ g ^ 1, the orbital period was even shorter. 
The current 5.57 hr orbital period might require an extended 
epoch of magnetic braking of the secondary (Rappaport et 
al. 1983; Taam 1983; Pylyser & Savonije 1988; and references 
therein). 

As discussed in the above references, magnetic brak- 
ing is expected only to bo effective for secondaries that do 
not have a fully convective envelope. Although strong mass 
transfer rates are possible before mass loss leads to such an 
envelope, this epoch is only expected to last lO^-lO'^ years 
(Rappaport et al. 1983). This is consistent with our previous 
estimates of the lifetime of the X 1822—371 system. We note 
that during this strong mass transfer phase, one generally 
expects the secondary to exceed its main sequence radius 
(Taam 1983; Kalogera & Webbink 1996). If we define the 
ratio of the secondary's Roche lobe radius to its main se- 
quence radius as 71, then given the discussion of §3 we can 
approximate the mass of the primary as 

M,^3.0n-'/'(—J^ll4^^]"'' . (10) 
\^70 km s-Vsin85°y ^ ^ 

Thus, 72. ^ 1.2 allows for a neutron star mass ^ 2 Mq. 

Finally we note that in addition to any mass loss asso- 
ciated with magnetic braking, further strong mass loss from 
such a system as suggested by King & Begelman (1999) leads 
to another explanation (in addition to postulating that the 
coronal optical depth t <^ 1 or t S> 1) for the fact that we 
observe only « 10~^ of the accretion luminosity inferred for 
an efficiency of r; ~ 10%. 'ADIOS' models wherein most of 
the accreted mass never reaches the surface or event horizon 
of the primary (Blandford & Begelman 1999; King & Ritter 
1999) allow for either a neutron star or black hole primary. 
In such a scenario, the energy of viscous dissipation is not ef- 
ficiently radiated away from the system and must be carried 
away via a wind. Viscous dissipation occurs out to the circu- 
larization radius, which is the inferred radius of the corona 
from our fits of §3, and thus dissipation may in part be re- 
sponsible for the presence of the corona (which in this case 
could be an ADIOS wind). 



6 CONCLUSIONS 

We have presented observations of the low mass X-ray bi- 
nary X 1822-371 taken with both RXTE and ASCA. We 
considered two broad band fits, which we took as approx- 
imately representing 'optically thick' and 'optically thin' 
coronal emission. Either model fit the data nearly equally 
well. Likewise, the X-ray lightcurves folded on the orbital pe- 
riod were also equally well-fit by optically thick or thin mod- 
els, and furthermore these fits could not distinguish among 
a white dwarf, neutron star, or black hole primary. High 
spectral resolution ASCA data revealed complex structure 
in the Fe Kq/K/3 region, consisting of possibly two lines. 
These latter features might be related to obscuring material 
between the edge of the corona and the disc rim, and also 
possibly related to emission from any optically thick regions 
of the corona. 

The spectral ambiguities of the line region will likely 
be resolved by upcoming high resolution observations with 
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XMM-Newton and Chandra. The nature of the primary, 
however, will likely not be revealed by these observations. 
As discussed above, understanding the nature of the pri- 
mary will in large part depend upon further observations 
of the secondary, and depend upon identifying a plausible 
evolutionary scenario for the X 1822—371 system. 
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APPENDIX A: DATA ANALYSIS 
METHODOLOGY 

Al RXTE Data Analysis 

We extracted data from both pointed instruments on RXTE, 
the Proportional Counter Array, PCA, and the High Energy 
X-ray Timing Experiment, HEXTE. X 1822-371, however, 
is both very faint and very soft, therefore we only spectrally 
fit data from the PCA instrument. The RXTE data were 
analyzed using the same procedure as that for our analysis 
of the spectrum of GX 339-4 (Wilms et al. 1999b). Specifi- 
cally, all RXTE results in this paper were obtained using the 
standard RXTE data analysis software, ftools version 4.2, 
and response matrix v3.1. Data selection criteria were that 
the source elevation was larger than 10° above the earth 
limb and data measured within 30 minutes of passages of 
the South Atlantic Anomaly or during times of high parti- 
cle background (as expressed by the "electron ratio" being 
greater than 0.1) were ignored. To increase the signal to 
noise level of the data, we restricted the analysis to the first 
anode layer of the proportional counter units (PCUs) where 
most source photons arc detected (the particle background 
is almost independent of the anode layer), and we combined 
the data from all five PCUs. We only used data where all 
five PCUs were turned on, which was nearly the entire ob- 
servation. 

For spectral fitting, we limited the energy range of 
the PCA data from 3 to 30keV To take into account the 
calibration uncertainty of the PCA we applied the chan- 
nel dependent systematic uncertainties described by Wilms 
et al. (1999b). These uncertainties were determined from 
a power-law fit to an observation of the Crab nebula and 
pulsar taking into account all anode chains; however, they 
do also provide a good estimate for the first anode layer 
only since most of the photons are detected in this layer. 
Background subtraction of the PCA data was performed 
using the 'SkyVLE' model, as for our previous studies of 
GX 339-4 (Wilms et al. 1999a). 



A2 ASCA Data Extraction 

We extracted data from the two solid state detectors {SISO, 
SISl) and the two gas detectors {GIS2, GIS3) onboard 
ASCA by using the standard ftools as described in the 
ASCA Data Reduction Guide (Day et al. 1998). We chose 
circular extraction regions with radii of fs 4arcmin for the 
SIS detectors, and ~ 6 arcmin for the CIS detectors. We 
excluded approximately the central 1 arcmin to avoid the 
possibility of photon pileup. We used the sisclean and 
gisclean tools (with default values) to remove hot and flick- 
ering pixels. We filtered the data with the strict cleaning 
criteria outlined by Brandt et al. (1996); however, we took 
the larger value of 7 GeV/c for the rigidity. We rebinned the 
spectral files so that each energy bin contained a minimum of 
20 photons. We retained SIS data in the 0.5 to 10 keV range 
and GIS data in the 1 to 10 keV range. The background was 
measured from rectangular regions on the two edges of the 
chip farthest from the source {SIS data), or from annuli with 
inner radii > 6 arcmin {GIS data). These data were cleaned 
and filtered in the same manner as the source files. 

Note that for the simultaneous ASCA/ RXTE obser- 
vations, we combined the two SIS detectors into a single 
spectrum, and we combined the two GIS detectors into a 
single spectrum, properly weighting the response matrices. 
We accounted for the cross-calibration uncertainties of the 
SIS and GIS instruments relative to each other and relative 
to RXTE by introducing a multiplicative constant for each 
detector in all of our fits. Note also that the phase filter- 
ing option in xselect produces flawed results, so the time 
filtering option was used instead to produce phase selected 
spectra. 
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